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Abstract: we give an approximate algorithm of computing holonomic sys- 
tems of linear differential equations for definite integrals with parameters. We 
show that this algorithm gives a correct answer in finite steps, but we have no 
general stopping condition. We apply the approximate method to find differ- 
ential equations for integrals associated to smooth Fano polytopes. They are 
interested in the study of K3 surfaces and the toric mirror symmetry. In this 
class of integrals, we can apply Stienstra's rank formula to our algorithm, which 
gives a stopping condition of the approximate algorithm. 

1 Introduction 

Let D be the ring of differential operators with polynomial coefficients. A func- 
tion /(xi, . . . , Xn) is called a holonomic function when the annihilating ideal in 
£) of / is a holonomic ideal (see, e.g., [ITj chapter I]). Finding systems of differ- 
ential equations for the definite integral F{xm+i, ■ ■ ■ , Xn) — f{x)dxi , . . . , Xm 
where C is a cycle is a fundamental problem in the symbolic computation. Since 
the celebrating work of D.Zeilberger about 20 years ago, a lot of algorithms have 
been proposed. Among them, the algorithm given by Oaku |13| constructs gen- 
erators of (/ -I- diD + • • • -f- dmD) n D' for a holonomic left ideal / in D. Here, 
D' = K{xrn+i^ ■ ■ ■,Xn,dm+i, ■ • • , 9„) and K is C. 

We give an approximate variation of this algorithm, which improves the 
performance of this algorithm. We show that the approximate algorithm gives 
a correct answer in finite steps, but we have no general stopping condition. 
We apply this method to find differential equations for integrals associated to 
smooth Fano polytopes. They are interested in the study of K3 surfaces and the 
toric mirror symmetry [l],[6], [7], [lOj . For this class of integrals, we can apply 
Stienstra's rank formula jl8j to our algorithm, which gives a stopping condition 
of the approximate algorithm. 
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2 Definite Integrals associated to Smooth Fano 
Polytopes 



A polytope V is called d-dimensional smooth Fano polytope when it satisfies the 
following 5 conditions. 

1. It is d-dimensional lattice polytope. 

2. The origin is in the interior of the polytope. 

3. The dual polytope V* is a lattice polytope. 

4. All faces are simplices (simplicial polytope) 

5. The vertices of each facet is a Z basis of Z"*. 

A polytope is called reflexive polytope when it satisfies the conditions 1, 2, 

3. 

A classification of reflexive polytopes in dimensions 3 and 4 is given in [7], 
[S] . The following list is the list of smooth Fano polytopes for dimensions 2 and 
3. It can be obtained by extracting smooth Fano polytopes from Kreuzer's list 
of reflexive Fano polytopes. We note that 0bro(|T51) gave an algorithm of find 
all smooth Fano polytopes and also gave a list of them up to the dimension 7 
modulo isomorphism. The next list is generated by the program by 0bro. 



dim 


index 


vertices 


2 





(1,0), (0,1), (-1,-1) 


2 


1 


(1,0), (0,1), (-1,0), (0,-1) 


2 


2 


(1,0), (0,1), (-1,1), (0,-1) 


2 


3 


(1,0), (0,1), (-1,1), (1,-1), (-1,0) 


2 


4 


(1,0), (0,1), (-1,1), (1,-1), (-1,0), (0,-1) 


3 





(1,0,0), (0,1,0), (0,0,1), (-1,-1,-1) 


3 


1 


(1,0,0), (0,1,0), (0,0,1), (-1,0,0), (0,-1,-1) 


3 


2 


(1,0,0), (0,1,0), (0,0,1), (-1,-1,1), (0,0,-1) 


3 


3 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,-1) 


3 


4 


(1,0,0), (0,1,0), (0,0,1), (-1,-1,2), (0,0,-1) 


3 


5 


(1,0,0), (0,1,0), (0,0,1), (-1,-1,2), (0,1,-1), (0,0,-1) 


3 


6 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (1,0,-1), (-1,-1,0) 


3 


7 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,0), (0,0,-1) 


3 


8 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,0,-1), (0,-1,-1) 


3 


9 


(1,0,0), (0,1,0), (0,0,1), (-1,0,0), (0,-1,0), (0,0,-1) 


3 


10 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,1), (0,0,-1) 


3 


11 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,1,-1), (0,-1,0) 


3 


12 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,1), (0,1,-1), (0,-1,0) 


3 


13 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,1), (0,1,-1), (0,0,-1) 


3 


14 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,1,-1), (0,-1,0), (0,0,-1) 


3 


15 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (1,0,-1), (-1,0,0), (0,-1,0) 


3 


16 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (0,-1,1), (0,1,-1), (0,-1,0), (0,0,-1) 


3 


17 


(1,0,0), (0,1,0), (0,0,1), (-1,0,1), (1,0,-1), (-1,0,0), (0,-1,0), (0,0,-1) 
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Let {ai, • • • ,a„i} be the vertices of a n-dimensional smooth Fane polytope 
Pn,k and Um+i the origin. The suffix k is the index in the hst of the smooth 



Fano polytopes. We put fn,k{x, t) = X]"=t^ ■ ^"^^ interested in a definite 
integral with parameters associated to these polytopes defined by 



Jc 



where a; is a generic parameter vector and C is a cycle in Hn{Tx,C), = 
{t G (C*)" I f{t,x) 0}. The definite integral is called a period integral 
associated to Pn,k- The function F^^k satisfies A-hypergeometric system for 
A = (ai,... , On+i) and /3 = (—1, 0, . . . , 0)"^ and the system is reducible [T]. The 
construction of the subsystem is implicit in [T] . We are interested in an explicit 
expression of the subsystem of the A-hypergeometric system satisfied by Fn^k- 

Example 1 The vertices of Pg^o are {(1, 0, 0), (0, f , 0), (0, 0, f ), (-f , -f , -f )}. 
Then, we have, by definition, 

Fzfl{x) = / [Xiti + X2t2 + X^ts + Xit^^t2^t^^)^^t^^t2^t^^dtidt2dtz, 

Jc 

which is nothing but the ^-hypergeometric integral for the matrix 

1 

1 n n - 



(I 1 1 

1 

1 

\0 1 



1\ 






and /3 = (— 1, 0, 0, 0)-^. It satisfies the ^-hypergeometric system Ha{P) of which 
holonomic rank is 4. As we will see later, the function F3fi{x) satisfies a system 
of which rank is 3. 

We note that the problem of finding systems of differential equations for 
period integrals associated to P„,fe has been considered to study the moduli space 
of the family of hypersurfaces f{t,x) = 0. Ishige |^ studied the case of P3 = 
{(1,0,0),(0, 1,0), (0,0, 1),(0,0,-I),(-I,-1,-I)} and Nagano [10] studied the 
cases of 

P2 = {(1,0,0),(0,1,0),(0,0,1),(-1,~1,0), (0,0,-1)}, 
P4 = {(1, 0, 0), (0, 1, 0), (0, 0, 1), (0, 0, -1), (-1, -1, -2)}, 
P5 = {(1,0,0), (0,1,0), (0,0,1), (-1,-1,0), (-1,-1,-1)}, 
P, = {(1,0,0), (0,1,0), (0,0,1), (0,-1,-1), (-1,0,-1)} 

The correspondence of our table and Nagano's table is as follows. 



P. :(dim,index) P2 : (3,1) P3 : (3,2) P4 : (3,4) P5 : (3,3) 



Pr is not simplicial polytope. Nagano gives series expansion of the function 
Fn^k and gives a heuristic method to find annihilating operators of the function 
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by the method of undetermined coefficients and this method is enough for his 
purpose. 

We wiU give an efficient algorithm to derive systems of differential equa- 
tions for period integrals in this paper. Our method gives a basis of associated 
cohomology groups. 



3 Computational Bottleneck of Integration Al- 
gorithm for L>-modules 

Computation of differential equations satisfied by a definite integral with pa- 
rameters can be reduced to the I?-module theoretic integral of the annihilating 
ideal of the integrand ([13], [S], [II])- This method works for any integrand 
which is a holonomic function. In particular, we can apply it to find differen- 
tial equations for the function Fn^k- However, it requires huge computational 
resources in general. In this section, we will explain what are bottlenecks of this 
algorithm in case of our problem. 

We denote by D the ring of differential operators of polynomial coefficients 
in n variables 

D = K{xi, . . . ,Xm,Xm+l,- . ■Xn,dl,.. . , dm,dm+l, • . • , 9„) 

and by D' that in n — to variables 

D K (^Xjy^^i , . . . , Xyi , drn+l : • ■ • : ^n) 

Let / be a left ideal of D. The integration ideal of / with respect to a;i, • • • , Xm 
is 

J= {I + diD + --- + d,nD)nD\ 

which is known to be a holonomic left ideal of D'. 

When the ideal / is the annihilating ideal of a function f{xi, . . . , Xn), the 
integration ideal J annihilates the definite integral 

f {xi, ' ' ' 7 Xmi 7 ^m+l 5 ' ' ' 5 Xji)dx\ • • ' dXm 

C 

for any m-cycle C (see, e.g., [TTl Chap5 ]). 

Let us review the Z?-module theoretic integration algorithm following Oaku 
[IB] to show bottlenecks of this algorithm. We define the ring isomorphism of 
DJFby 

I Xi (to < I <n) \Oi [m < I < n) 

and we call it the Fourier transformation. The inverse map of F is called the in- 
verse Fourier transformation and denoted by . We use the symbol \'a.(^-w,w) 
defined in |I7| to express the initial ideal. The generic 6-function of / with re- 
spect to the weight w G R" is the monic generator of the ideal in(_^_u,)(/)nC[s] 
in C[s] where s = Yll=i WiXidi. 
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Algorithm 1 (Computation of integration ideal, |13| . |14] . |17| ) 

Input: A set of generators of holonomic left D ideal I . 

A weight vector w — (wi, . . . , Wm, u'm+i, ■ ■ ■ ,Wn) satisfying 
,...,«;,„> 0, w„i+i ^ ■■■ = Wn = 0. 

Output: A set of generators of the integration ideal of I with respect to Xi, . . . , Xm- 

1. Compute the restriction ideal of with respect to the weight vector w. 
This step consists of the following steps. 

(a) Compute a Crohner basis G = {hi, ■ ■ ■ ,hi} of with respect to 

the order <(-w.w)- 
(h) Compute the generic h-function b{s) of with respect to {—w,w). 

If there is no non-negative integral root of b{s) = 0, then output 1. 

Let sq be the maximal non-negative integral root of b{s) = 0. 

(c) Put nil ^ ords^_w,^){hi), 

Bd = {di • • • 9^ I iiwi H h imWm < d}, 

^ = • • ■ , «m) I hWl H h imWm < Sq} . 

(d) PutB = {(a^/i.)Ui=...=.„=o I e Bs,-,n^,h, G G}. Here, (a'5/i,)Ui=-=.™=o 
is calculated by substituting Xi — ■ ■ ■ — Xm = after expanding d^hi 

into the form X]« u '^".-"^"^"('^".i' ^ constant). 

2. Elements in J-^^{B) can be expressed as 

wiii H \-Wjni^ri 

We regard it as an element of (D'Y . In other words, we consider a vector 

(cij.... ,i„)toiiiH hiumim<(i '^f which entries belong to D' . Let M be the left 

D' module generated by them. 

3. Compute a Crobner basis of M with respect to an POT order such that the 
element standing for x^ - ■ ■ cc^ is the minimum. Collect elements in the 
Crobner basis such that all entries except the entry standing for xi ■ ■ ■ 

is 0. We regard these elements as elements of D' . Output them. 

Example 2 (Period integral associated to a 2 dimensional smooth Fano polytope) 

We will derive a system of differential equations for the period integral 

F2.o{x) — / {xiti + X2t2 + X3ti^t2^ + Xi)^^ti^t2^dtldt2 

Jc 

associated to the smooth Fano polytope P2.0 algorithm we have presented. 

1. Obtain the annihilating ideal I of 

f2,o{x, t) = {xiti + X2t2 + xstih^^yHih^^ 

by Oaku's algorithm (1121). 
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2. For the ideal I, we compute the integration ideal of I with respect to ti , t2 



J ^ {I + dt,D + dt^D)nD' 

Here, D = K{xi,X2,x^,Xi,ti,t2,dx^,dx2,dx3,dxi,dt^,dt^) , 
D' = K{xi,X2,X'i,Xi,dx^,dx2,dx3,dxi). 

Generators of J obtained by Algorithm]^ are 

{xl + 27xiX2Xz)d'i^ + ixldxi + X4, 9x2X^81^ — x\dxidxi - x^dxi, ^xix^di^ — x\dx2dxi — XAdx2, 

2 2 2 2 2 

-%xiX2dx^+Xidx3dxi+X4,dx3, -3x3dx^ ~ X4dxidx2, ~ 'ix2dx^ ~ Xidx^dx^, ~ i^xidx^ - Xidx^dxs, 
- di^ + dxidx2dx3, Xidxi + ixxdx^ + 1, - xaOx^ - Sx2dx2 - 1, Xidx^ + Sx^dx^ + 1 

The four elements in the last line are elements in the A-hypergeometric ideal 

/I 1 1 1\ 

Ha{,P) and other elements are not. Here, we put A ^ \1 —1 ,/3 = 

\0 1 -1 0/ 

(^i,or. 

We present timing data obtained by applying these algorithms to our prob- 
lem. The timing data are taken on Risa/Asir (111]) version 20100206 and the 
library package nkjrestriction.rr running on a machine with Xeon X5470 
(3.33GHz) CPU and 3.6 GB memory except the computation of annihilating 
ideal, which is computed on Singular and the same machine. It is known that 
the Singular implementation of it is the best one among several implementa- 
tions. In the table below, the field Ann stands for Briangn-Maisonobe algo- 
rithm of computing annihilating ideal ([5], [H]) implemented in Singular, the 
field generic-b for the steps 1 (a), (b) of Algorithm [U the field base for the steps 
1 (c) , (d) and the step 2 of Algorithm [TJ and the field GB for the step 3 of 
Algorithm[T] Each entry means seconds. The symbol " — " means that we could 
not obtain the result in one day. 
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dim 


index 


Ann 


j fjj yj ] V.J 1 J CAiJ.i-V-1 ^VjJ.i.v-i. u 




GB 


2 







0.004 


0.04 


0.004 


2 


1 




0.022 


0.04 


0.008 


2 


2 




0.035 


0.03 


0.019 


2 


3 


2 


0.58 


0.11 


0.19 


2 


4 


180 


93 


1.6 


14 


O 


n 


I 


052 


048 


020 


3 


1 


1 


0.20 


0.09 


0.052 


3 


2 


3 


0.36 


0.16 


0.10 


3 


3 


5 


1.6 


0.25 


0.14 


3 


4 


14 


1.4 


0.25 


0.013 


3 


5 


4165 


1710 


726 


5726 


3 


6 


2383 


5037 


1007 


6194 


3 


7 


135 


92 


87 


257 


3 


8 


1551 


1618 


499 


2227 


3 


9 


16 


13 


34 


86 


3 


10 


1492 


852 


183 


994 


3 


11 


6213 


1588 


406 


1839 


3 


12 











The data tell us that if we can accelerate the steps "Ann" (computation of 
annihilating ideal) and "<(_t„^^,)-GB and generic-b" (computation of {—w,w)- 
Grobner basis and generic 6- function), then we may be able to speed up the 
computation in this list. In the next section, we will suggest a new efficient 
method. 

4 Approximate Integration Algorithm 

Definition 1 (Approximate annihilating ideal) 

Let /, g be polynomials. The i-th approximate annihilating ideal of the rational 
function ^ is the ideal generated by the elements of which (0, l)-order is less 

than or equal to i in Ann/j ^ . We denote it by Ann^'' ^ . 
The following method is simple, but is very useful. 

Algorithm 2 (Computation of approximate annihilating ideal [3]) 

Input: rational function ^, order i 

Output : generators of i-th approximate annihilating ideal Ann^'' ^ . 

1. Put P — aad" where Oq, are polynomials to be determined and the sum 
runs over all a such that the order of {ai + • • • + Q!„ < i) is less than 
or equal to i. 
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2. Apply the differential operator P to the rational function |- and derive a 

relation for polynomials Ua so that P' ^ — 0. The polynomials can be 
determined by a syzygy computation in the ring of polynomials. 

Example 3 (Computation of approximate annihilating ideal) 

We will compute the first approximate annihilating ideal for 



1. Apply P = aidti + a2dt2 + a^d^^ + a^dx^ + a^dx^ + a^dx^ + 07 (a^ are 
polynomials to be determined) to f2.o{tix)- Then, the numerator n{t,x) 
is 

n{t,x) ={—t2X4 — t\x2 — 2tit2Xi)ai + {—tix^ — 2tit2X2 — t\xi)a2 — 

t\t2az - tit\ai - as - iit2a6 + {tit2Xi + 2:3 + titlx2 + tjt2Xi)a7 

2. We want to determine so that n{t,x) — 0. Let Ci be the coefficient 
of Qi of n{t, x). Any element of the syzygy module Syz(ci, • • • ,07) gives 
{ai} such that n(t,x) — and if a set {a^} makes n(t,x) — 0, then it is 
an element of the syzygy. The syzygy can be obtained by Grobner basis 
computation in the ring of polynomials and it is generated by 

(-1, 0, 0, 0, ^22:4 + t22:2, 2a;i, 0), (0, -1, 0, 0, ^1x4 + i^xi, 2x2, 0), (0, 0, -1, 0, 0, h, 0), 
(0, 0, 0, -1, 0, t2, 0), (0, 0, 0, 0, ~tit2, 1,0), (0, 0, 0, 0, X3, X4 + t2X2 + hxi, 1), 
(0,0, 0,0, 0,tit2a;4 + X3 + tit\x2 + t\t2Xi,tit2) 

The set of differential operators standing for them 

- 9*1 + (t2X4: + tlx2)dx^ + 2xidxi,-dt^ + (tia;4 + tlxi)dx3 + 2x28^^,-8^;^ + hd^ 

- + t2dxi,-tit2dxs + dx^.Xsdxs + {x^ + t2X2 + tiXi)dx4 + 1, 
(tit2X4 +Xz+ tit\x2 + t\t2Xi)dxi + tit2 

is a set of generators of the first approximate annihilating ideal. 

It is an interesting open question that at which order the approximate an- 
nihilating ideal equals to the annihilating ideal. Castro and Ucha [5] proved 
that if a polynomial / in two variables is weighted homogeneous, then we have 
Ann^^ J — Ann/3 J. 

Utilizing approximate annihilating ideal, we give an approximate integration 
algorithm. 

The syzygy computation in the ring of polynomials is faster than computa- 
tion of the annihilating ideal in D in general. Then, we can expect that our 
approximate integration algorithm runs faster than the standard one. However, 
we need a stopping condition to stop the approximate procedure. 
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Algorithm 3 (Approximate integration algorithm) 

Input : ^ is a rational function, m is a natural number. 

Output : A subideal J' of the integration ideal of Ann^i^. 

1. Compute J = Ann^"^i (Algorithm [2]) 

2. If J is holonomic, then apply the (Z)-module) theoretic integration algo- 
rithm to J and put it J', else put J' = (0). 

Theorem 1 There exists a natural number m such that the approximate inte- 
gration ideal J' agrees with the integration ideal. 

Proof. Since D is a Noetherian ring, there exists m such that Ann^*' ^ = 
Anno-. The theorem follows from this fact. Q.E.D. 



Example 4 (m is 2.) 

We use the example by Castro and Ucha([3]). Consider the polynomial 
f = + + xy'^ (Reiffen curves). We have 

Ann^^^ J C Anng^y = Ann^y 

The first approximate integration ideal for j with respect to x is J^^'^ = 
D ■ {yP} and the integration ideal J is 13 • {P}. Here, 

P = (-27/+256y3)9,y + (-432y3+3456y2)a2+(-18962/2+12336y)aj,-2184y+10920 
For m — 2, the approximate integration ideal J^^^ agrees with J . 

Let us come back to our problem of computing differential equations for pe- 
riod maps associated to smooth Fano polytopes. In application, we do not need 
to find the exact integration ideal and we may find an approximate integration 
ideal of which holonomic rank agrees with that of the exact solution. For this 
purpose, the method given by Stienstra gives a way to find the approximate 
order "m" . We follow the notation of the paper of Stienstra [181 P-434] 

T\Zs is defined as C""^\{a;i • • • Xn-vf{xi, ■ ■ ■ ,a;„_i) = 0} where f {xi, ■ ■ ■ ,Xn-i) 

is defined as follows from the set of points A—\^ ^ | . 

\ai 02 ■■■ aij 

I 

/(xi, • • • ,a;„_i) = ^Uix"* (ui e C) 

i=l 

Algorithm 4 (An algorithm finding a set of differential operators for integrals 
associated to a class of smooth Fano polytopes) 
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1. Evaluate the dimension r — dhxiW„H'^ '^(T\ Z^) by Stienstra's method 
(UM p435 (57), p448 Th.lO (iv)]). 

2. i = l. 

3. Compute the integral ideal J for Ann^'' (/) . 

^. If the holonomic rank J is equal to r, then stop and output J, else i i + 1 
and goto the step 3. 

Theorem 2 

1. Suppose that A is smooth Fano (reflexive Gorenstein) and admits a uni- 
modular regular triangulation and, moreover. 



Proof. (1) Since the annihilating ideal J — Aimjj (/) is finitely generated, 
there exists i such that D ■ Ann^''(/) = Ann^i (/). The integration module of 
D/J is isomorphic to if"~^(T \ Z^) [M]. Therefore, our algorithm outputs the 
integration ideal of which rank is equal to r. (2) We have W^„„iif "^^(T\ Z^) C 
W„i?"-i(T \ Zs) C F"-i(T \ Zs) in general. The claim (2) follows from this 
inclusion. Q.E.D. 



Example 5 We consider the case P2.o- Vertices of the polytope P2.0 are ai = 
(1,0), 02 = (0,1), 03 = (-1,-1). We put 04 = (0,0) and 



Put r = {1, 2, 3, 4, 12, 13, 23, 124, 234, 134}. Here, ii---is denotes the simplex 
spanned by a^^ , • • • , a^^ . Then, 7~ is a unimodular triangulation of the set of 
points {ai,...,a4}. Following Stienstra |18| . we consider the ring TZa,t = 
Z[ci, C2, C3, C4]/j7 where Ci is a variable standing for the point Oi and the ideal 
J' of the ring Z[ci, C2, C3, C4] is 



The first 3 linear polynomial stand for A and the last monomial stands for a sim- 
plex which is not in the triangulation T. Put coreT — {intersection of the maximal simplices} 
and Ccore = Y\i£coreT ^^^^ example, Ccore — C4. Stienstra's theorem claims 

that the rank of TZa,t / ^^^tZa r'^core as Z-module gives a lower bound of the 
holonomic rank of the approximate integration ideal. The first approximate in- 
tegration ideal gives this lower bound for this example and then we obtain the 
answer. 



r = diniW„iJ""i(T\ Z,) = dim iJ"-i(T \ Z,) (1) 

Then, Algorithm outputs a set of differential operators of which holo- 
nomic rank agrees with r. 

2. If Algorithm^ stops, then it implies the equality (1^. 




J = {C1+C2+C3+ C4, Ci - C3, C2 - C3, C1C2C3) 



10 



Wc have applied this algorithm and obtain the following two tables. The field 
"tank of HA{Py^ is the holonomic rank of the hypergeometric system Ha{(3) 
and the field "rank of W„i?"~^" is the lower bound evaluated by Stienstra's 
method. The field "AppAnn" in the second table is the time to find the first 
order approximate annihilating ideal. It is an interesting observation that we 
need only first order operators. 



dim 


index 


rank of Ha{I3) 


rank of W„i/""^ 


2 





3 


2 


2 


1 


4 


2 


2 


2 


4 


2 


2 


3 


5 


2 


2 


4 


6 


2 


3 





4 


3 


3 


1 


6 


4 


3 


2 


6 


4 


3 


3 


6 


4 


3 


4 


6 


4 


3 


5 


8 


5 


3 


6 


8 


5 


3 


7 


8 


5 


3 


8 


8 


5 


3 


9 


8 


5 


3 


10 


8 


5 


3 


11 


8 


5 


3 


12 


10 


6 


3 


13 


10 


6 


3 


14 


10 
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5 Conclusion 



We have succeeded to efficiently find differential equations for period maps 
standing for P3.5 to P^^u by our new approximate method. Explicit expres- 
sions of differential equations for period maps and programs used in this paper 
can be found in |http : //www . math . kobe-u . ac . jp70 penXM/Math/sf eLno| 

Let us briefly discuss two future research plans of this work. 

The first one is to utilize the construction of all irreducible quotients of A- 
hypergeometric D-module by Saito [1^. To do this, we need to give a new 
integration algorithm to fit to his construction. 

The second one is to improve the final step of computing a <(_i„^iL,)-Grobner 
basis. Let R be the ring of differential operators with rational function coeffi- 
cients with the variables xi, . . . , Xm, Xm+i, . . . , a;„ and R' the ring of differential 
operators with the variables Xm+i, ■ ■ ■ ,Xn- Chyzak A\ gives an approximate 
algorithm to find elements of {RI + X]t=i diR) n R' by utilizing the method 
of undermined coefficients. Since we have to find an approximate basis of the 
cohomology group to use Stienstra's criterion, we cannot use the R and we need 
to use D. By this reason, we cannot try this idea with the current implemen- 
tation mgf un, which uses R, of this algorithm in the Maple. However, Chyzak's 
algorithm itself can be easily modified for D. It is our future project to give an 
efficient implementation of this approximate algorithm for D to apply to our 
problem. 
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